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We summarize our recent results on the phase diagram of QCD with Nf = 2-\-l quark flavors, as 
a function of temperature T and quark chemical potential jj.. Using staggered fermions, lattices 
with temporal extent Ni ~ 4, and the exact RHMC algorithm, we first determine the critical line in 
the quark mass plane (m„ ;/,ms) where the finite temperature transition at = is second order. 
We confirm that the physical point lies on the crossover side of this line. Our data are consistent 
with a tricritical point at {m^j^nis) = (0, ^ 500) MeV. 

Then, using an imaginary chemical potential, we determine in which direction this second-order 
line moves as the chemical potential is turned on. Contrary to standard expectations, we find 
that the region of first-order transitions shrinks in the presence of a chemical potential, which is 
inconsistent with the presence of a QCD critical point at small chemical potential. 
The emphasis is put on clarifying the translation of our results from lattice to physical units, and 
on discussing the apparent contradiction of our findings with earlier lattice studies. 
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Figure 1: Schematic phase transition behaviour of A'y = 2 + 1 flavor QCD for different choices of quark 
masses (m„ f/,ms), at = (from [|l|]). 

1. Introduction 

In recent years, considerable efforts have been devoted to the determination of the phase di- 
agram of QCD at finite temperature and density [|l|]. At zero chemical potential, the behaviour as 
a function of temperature depends on the quark masses mj, ^ and wis, and expectations are summa- 
rized in Fig. 1. In the limits of zero and infinite quark masses (lower left and upper right corners), 
order parameters corresponding to the breaking of a symmetry can be defined, and one finds nu- 
merically that a first-order transition takes place at a finite temperature Tc. On the other hand, 
one observes an analytic crossover at intermediate quark masses. Hence, each corner must be sur- 
rounded by a region of first-order transition, bounded by a second-order line as in Fig. 1. The line 
in the heavy-quark comer has been studied in [^]. Here, we want to determine the chiral critical 
line. 

Along both lines, the universality class is that of the 'id Ising model. Therefore, a powerful 
tool to determine the critical couplings is the Binder cumulant B4 = {SX"^) / {SX^Y, where 5X = 
X — {X), and we take for X the u,d quark condensate yy- On lattices 8^, 12^ and 16^ x 4, we 
estimate the critical couplings as those for which B4 = 1.604.., the 3d Ising value. For each mass 
point {mu^d,ms), we accumulate at least 200k RHMC trajectories, and interpolate among 4 or more 
rriiij values to find the critical rriuj mass m'^ at a given m^. We obtain the set of points in Fig. 4, left. 

We then consider the effect of a baryonic chemical potential. As a function of /i, represented 
vertically in Fig. 2, the critical line determined at /i = spans a surface. The standard expectation 
is depicted in Fig. 2 left. The first order region expands as n is turned on, so that the physical point, 
initially in the crossover region, eventually belongs to the critical surface. At that chemical potential 
He, the transition is second order: that is the QCD critical point. Increasing /i further makes the 
transition first order. A completely different scenario arises if instead the first-order region shrinks 
as /I is turned on. In that case (Fig. 2 right), the physical point remains in the crossover region for 
any jj.. Since the phenomenologically interesting question is whether a QCD critical point (He^Te) 
exists at small pletP^e/Te ^ 1, this question can be addressed by an analytic continuation based 
on a Taylor expansion. Using an imaginary chemical potential [||, ^, we determine the curvature 
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Figure 2: The chiral critical surface in the case of positive (left) and negative (right) curvature. If the 
physical point is in the crossover region for /i = 0, a finite /i phase transition will only arise in the scenario 
(left) with positive curvature, where the first-order region expands with Note that for heavy quarks, the 
first-order region shrinks with |^| (right) [|]]. 



^1^=0 of the critical surface at /i = 0. We find that it is negative, so that the first-order region 
shrinks as in Fig. 2 right. Note that in the opposite corner, the first-order region also shrinks 

Sec. 2 tests our methodology in the Nf = 3 case. Sec. 3 describes the Nf = 2 + I study. 
Sec. 4 compares our results with earlier lattice studies and discusses the various limitations of our 
approach. 



2. A^^ = 3 

We first check our methodology in the case of 3 degenerate flavors. This is basically a repeat 
of Ref. this time using the RHMC algorithm [0] instead of the R algorithm. 

The RHMC algorithm eliminates the stepsize error of the R algorithm, which differs in mag- 
nitude in the chirally symmetric and broken phases As a result, the value of m''(/i = 0) is 
considerably different: = 0)) moves from 0.033(1) (R alg.) [|] to 0.0260(5) (see Fig. 3, 

left). We have checked, by performing zero-temperature simulations at this quark mass, that this 
is not a simple renormalization effect, but that the physical ratio mji/Tc is lowered by about 10%. 
Therefore, an exact algorithm appears mandatory for the study of the Nf = 2 + 1 critical line. More- 
over, RHMC turns out to be vastly more efficient, by up to a factor 20 in our case for the smallest 
quark masses [^]. 

We now turn on an imaginary chemical potential n = ipLj, and for each /i/ monitor the Binder 
cumulant S4 as a function of the quark mass. Our results are summarized in Fig. 3, left. The 
chemical potential has almost no influence on B4. A lowest-order fit, linear in am and (ajj.)-^, gives 
the eiTor band Fig. 3, right, corresponding to 

am' (an) =0.0270(5) -0.0024{m){a^f (2.1) 

Care must be taken for the conversion to physical units. The crucial point is that, as we 
increase the chemical potential Hj, we tune the gauge coupling j8 upwards to maintain criticality. 
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Figure 3: Left: B4{am,apLi) for different imaginary chemical potentials. Right: One-sigma error band for 
the critical mass am'^{ajj.i) resulting from a linear fit. 



so that a(j3) decreases: our observation that am'^{lJ,[) ^ const, does not mean that nf{}Xj) const., 
but that m'^ {iJ.i) increases with /i/, or decreases with a real chemical potential /i. If we express 

am'^ (0) am'^ (0) 



then ci and Cj are related by 



where m = m'^{jJL) in the second term. Writing the transition temperature as 



Tc(m,Ll) m — n£, / U \^ 

= l+A —^+b{-^] +... (2.5) 



one obtains 

c^ = (b+'-^—\—\\\-A'^\ . (2.6) 

c\ and ^ are both small, so that c\ is nearly equal to B. Estimates of B and A can be obtained by 
converting our result for the pseudo-critical gauge coupling 

j8o(am,aAt) =5.1369(3) + 1.94(3)(am-am[;) + 0.781(7)(aAt)^ (2.7) 
to physical units. Using the 2-loop /3 -function gives A = 2. Ill (17), B = —0.667(6) so that finally 

The error above is conservative and includes the uncertainty from using different fitting forms (see 
Table 2, Ref. [10]). The main source of systematic error comes from using the 2-loop j8-function 
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Figure 4: Left: The chiral critical line in the bare quark mass plane at = 0. Nf — 3 \& shown by the 
solid line. Also shown are the physical point according to JI^], and a fit corresponding to a tricritical point 
m'™ ~ 2. sr. Right: Comparison of the critical line at /i = and apLi ~ 0.2. 



to obtain B. The non-perturbative j3 -function varies more steeply and may increase A and B, in 
magnitude, by up to a factor 2. This will make ci more negative. 

We thus have clear evidence that, in the Nf = 3 theory on anNt =4 lattice, the region of first- 
order transitions shrinks as a baryon chemical potential is turned on, and the "exotic scenario" of 
Fig. 2, right, is the correct one. This result is further supported by recent simulations of the same 
theory, under an isospin chemical potential [11]. 



3. A^^ = 2 + 1 

We now proceed to the non-degenerate case. First, at = 0, we map out the line of second- 
order transitions in the (am„ £/,ani.s) plane. Our results, shown Fig. 4, left, are in qualitative agree- 
ment with expectations Fig. 1. In particular, they are consistent with the possible existence of a 
tricritical point (m„ j = 0,ms = m'"'^). Using its known, Gaussian exponents, our data favor (blue 
line in Fig. 4 left) a heavy m'/"^ ~ 2.8rc. 

A more immediate issue is whether the QCD physical point lies on the crossover side of the 
critical line as expected. For that purpose, we have performed spectrum calculations at T ~ 0, at 
the parameters corresponding to the horizontal arrow in Fig. 4 left {aniu^d = 0. 005, am^ = 0.25, jS = 
5.1857). They show that is approximately tuned to its physical value ~ ^|phys)> while the 
pion is lighter than in QCD = 0.148(2) < 0.18). This confirms that the physical point lies on 
the right of the critical line, i.e. in the crossover region ^ This conclusion has been confirmed by 



very recent calculations on finer lattices [ 13]. Also, we find Tc to vary little along the critical line. 



in accordance with model calculations [14]. 



We now couple an imaginary chemical potential ajj./ = 0.2 to the two light flavors, and measure 
the change in the critical mass am^ j as in the Nf = 3 case. Fig. 4 right shows the same trend as 
for Nf = 3: the critical mass is constant or slightly increasing, in lattice units. The conversion to 



'in fact, our estimate of the lattice parameters corresponding to the physical point is consistent with that of Fodor 
& Katz using the same action, but the R algorithm [[l2||. 
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physical units proceeds as in eqs.( p^p^ ). Since the critical gauge coupling po^aHi) increases 
with Hj, the coefficient B, which is the dominant contribution to c\, is negative. Together with a 
very small or slightly negative value for c\ , it implies again that the first-order region shrinks as the 
baryon chemical potential is turned on, and the "exotic scenario" of Fig. 2, right, is the correct one. 

This statement comes with several caveats: (/) our lattice is very coarse (a ~ 0.3 fm); (//) 
as we consider lighter our box becomes small (nijiL ~ 1.7 for the worst case); (///) we use 
"rooting" of the staggered determinant to simulate 1 and 2 flavors, albeit our measure is positive 
with an imaginary n, so that we avoid the pitfalls of [ 15]. 



4. Discussion 

Our results appear in qualitative contradiction with those of Fodor & Katz [|l^ and of Gavai 
& Gupta [|l6|], which both conclude for the existence of a critical point {iJ-e,Te) at small chemical 
potential He/Te < 1. Let us consider the reasons for such disagreement. 




Figure 5: Left: Effect of keeping the quark mass fixed in lattice units in [[l2|]. Right: Comparison at finite jj. 
between the Nf = 2+1 and the Nf = 2 theory considered in JI^ ] . 

• Fodor & Katz obtain Monte Carlo results at j8 = jSc, /X = 0, and perform a double reweighting 
in (PjH) along the pseudo-critical line ^c{o-lx). By construction, this reweighting is performed at a 
quark mass fixed in lattice units: am^ ^ = -r^ = const.. Since the critical temperature Tc decreases 
as they turn on /x, so does their quark mass. This decrease of the quark mass pushes the transition 
towards first order, which might be the reason why they find a critical point at small pL. This effect 
is illustrated in the sketch Fig. 5, left, where the bent trajectory a la Fodor & Katz intersects the 
critical surface, while the vertical line of constant physics does not. 

Put another way, Fodor & Katz measure the analogue of eq.(2J) instead of ( |2.3[ ). From their 
Fig. 1 (Ref. 1 12]), the coefficient c\ which one would extract would be essentially zero like ours. 



As in our case, the variation of with }X makes a dominant contribution, which may change the 
results qualitatively. 

• Gavai & Gupta try to infer the location of the critical point by estimating the radius of 
convergence of the Taylor expansion of the free energy in {jx/nT)-^. Regardless of the systematic 
error attached to such estimate when only 4 Taylor coefficients are available, we want to point out 
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that they consider a theory without strange quark, i.e. Nf = 2 only. The {lJ.,T) phase diagram 
of such a theory is qualitatively different from that of A^y = 2 + 1 QCD. At = 0, the order of 
the finite-temperature transition as m„ rf ^ is not settled [[l7|]. Assuming a second-order 0(4) 
transition, one expects then a tricritical point at {muJ = = jJ."'"^), beyond which a non-zero 
critical mass m'^^{iJ.) can be defined, as sketched in Fig. 5 right. The quantitative relevance of 
results, even accurate, for this Nf = 2 theory to QCD is unclear to us. 

Therefore, we find no inconsistency between our results and those above. We conclude that 
the existence of a critical point {ij.e,Te) in QCD at small chemical potential /^fi/Ti; < 1 is an open 
question. Our numerical evidence, with the caveats mentioned in Sec. ^ is that the curvature of 
the critical surface is as illustrated Fig. 2 right. Our main systematic error comes from our coarse 



lattice spacing a ~ 0.3 fm [18]. If confirmed on a finer lattice, the implications of our finding 
are as follows. In the region where a leading Taylor expansion of the critical surface is a good 
approximation, i.e. jx/T < 1, corresponding to the experimentally accessible regime, no critical 
point exists which is analytically connected to /i = 0. Of course, we cannot exclude that the QCD 
phase diagram is more complex, and partly inaccessible to our imaginary jx + Taylor expansion 
strategy. 
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